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Abstract 

We present in detail the calculation of the 0{as) virtual corrections to the matrix element 
for b s'j. Besides the one-loop virtual corrections of the electromagnetic and color dipole 
operators O7 and Og, we include the important two- loop contribution of the four- Fermi operator 
O2. By applying the Mellin-Barnes representation to certain internal propagators, the result 
of the two- loop diagrams is obtained analytically as an expansion in mc/mb. These results are 
then combined with existing 0(0;^) Bremsstrahlung corrections in order to obtain the inclusive 
rate for B Xg'y. The new contributions drastically reduce the large renormalization scale 
dependence of the leading logarithmic result. Thus a very precise Standard Model prediction 
for this inclusive process will become possible once also the corrections to the Wilson coefficients 
are available. 
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1 Introduction 



In the Standard model (SM), flavor- changing neutral currents only arise at the one-loop level. 
This is why the corresponding rare B meson decays are particularly sensitive to "new physics" . 
However, even within the Standard model framework, one can use them to constrain the 
Cabibbo-Kobayashi-Maskawa matrix elements which involve the top-quark. For both these 
reasons, precise experimental and theoretical work on these decays is required. 

In 1993, B —>■ K*^ was the first rare B decay mode measured by the CLEO collaboration 
Recently, also the first measurement of the inclusive photon energy spectrum and the 
branching ratio in the decay B ^ Xs + '~f was reported 0. In contrast to the exclusive 
channels, the inclusive mode allows a less model-dependent comparison with theory, because 
no specific bound state model is needed for the final state. This opens the road to a rigorous 
comparison with theory. 

The data agrees with the SM-based theoretical computations presented in [|, ^, |[, given 
that there are large uncertainties in both the experimental and the theoretical results. In 
particular, the measured branching ratio BR{B Xg'j) = (2.32 ± 0.67) x 10~^ overlaps 
with the SM-based estimates in [|, || and in @, 0] . 

In view of the expected increase in the experimental precision, the calculations must be 
refined correspondingly in order to allow quantitative statements about new physics or standard 
model parameters. So far, only the leading logarithmic corrections have been worked out 
systematically. In this paper we evaluate an important class of next order corrections, which 
we will describe in detail below 0. 

We start within the usual framework of an effective theory with five quarks, obtained by 
integrating out the heavier degrees of freedom which in the standard model are the top quark 
and the VF-boson. The effective Hamiltonian includes a complete set of dimension-6 operators 
relevant for the process 6 — > 57 (and b — » sjg) |^ 



V2 j=i 



with Gf being the Fermi coupling constant and Cj(/i) being the Wilson coefficients evaluated 



at the scale fi, and Xt 
as follows: 



VtbV-i-l with Vij being the CKM matrix elements. The operators Oj are 



Oi = 




02 = 


{cLaYbLa) ( 




03 = 


{sLal^bLa) 




04 = 


(sLal^bils) 




05 = 


{sLal^bLa) 


{urp'^^.Urp) + . 


Oe = 


{sLal^bLp) 




Or = 


(e/167r2)s« a^"' {mb{fi)R + i 


Os = 




(T^'^ (mb(/i)i? + 



:i.2) 



^Some of the diagrams were calculated by Scares 
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In the dipole type operators O7 and Og, e and F^^ {g^ and G^^) denote the electromagnetic 
(strong) couphng constant and field strength tensor, respectively. L = (1 — 75)/2 and R = 
(1 + 75)72 stand for the left and right-handed projection operators. It should be stressed in 
this context that the explicit mass factors in O7 and Og are the running quark masses. 

QCD corrections to the decay rate for b —>■ s'j bring in large logarithms of the form 
a^{mw) ^og"^{mb/M), where M = rrit ot m\Y and m < n (with n = 0,1, 2, ...). One can system- 
atically resum these large terms by renormalization group techniques. Usually, one matches the 
full standard model theory with the effective theory at the scale mw At this scale, the large 
logarithms generated by matrix elements in the effective theory are the same ones as in the full 
theory. Consequently, the Wilson coefficients only contain small QCD corrections. Using the 
renormalization group equation, the Wilson coefficients are then calculated at the scale /i ~ mi,, 
the relevant scale for a B meson decay. At this scale the large logarithms are contained in the 
Wilson coefficients while the matrix elements of the operators are free of them. 

As noted, so far the decay rate for 6 — > 57 has been systematically calculated only to 
leading logarithmic accuracy i.e., m = n. To this precision it is consistent to perform the 
'matching' of the effective and full theory without taking into account QCD-corrections 



and to calculate the anomalous dimension matrix to order [|TT|. The corresponding leading 
logarithmic Wilson coefficients are given explicitly in |T^. Their numerical values in the 
naive dimensional scheme (NDR) are listed in table 1 for different values of the renormalization 
scale II. The leading logarithmic contribution to the decay matrix element is then obtained by 
calculating the tree-level matrix element of the operator C7O7 and the one-loop matrix elements 
of the four-Fermi operators CiOi {i = 1, 6). In the NDR scheme the latter can be absorbed 
into a redefinition of C7 Cj^^ 

Cff = CT + QdC, + ?,QdC, . (1.3) 

In the 't Hooft-Veltman scheme (HV) [1^ , the contribution of the four- Fermi operators vanishes. 
The Wilson coefficients C^ and Cg in the HV scheme are identical to Cj^^ and C^^^ in the NDR 
scheme. Consequently, the complete leading logarithmic result for the decay amplitude b ^ S'j 
is indeed scheme independent. 

Since the first order calculations have large scale uncertainties, it is important to take into 
account the next-to-leading order corrections. They are most prominent in the photon energy 
spectrum. While it is a (5-function (which is smeared out by the Fermi motion of the 6-quark 
inside the B meson) in the leading order, Bremsstrahlung corrections, i.e. the process b s'yg, 
broaden the shape of the spectrum substantially. Therefore, these important corrections have 
been taken into account for the contributions of the operators O7 and O2 some time ago 0] 
and recently also of the full operator basis |Q, |T1|, 0. As expected, the contributions of O7 
and O2 are by far the most important ones, especially in the experimentally accessible part 
of the spectrum. Also those (next-to-leading) corrections, which are necessary to cancel the 
infrared (and collinear) singularities of the Bremsstrahlung diagrams were included. These are 
the virtual gluon corrections to the contribution of the operator O7 for b sj and the virtual 
photon corrections to Og for b ^ sg. 

A complete next-to-leading calculation implies two classes of improvements: First, the Wil- 
son coefficients to next-leading order at the scale /i ~ ruh are required. To this end the matching 
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For the analogous b ^ sg transition, the effects of the four-Fermi operators can be absorbed by the shift 

Cs->CI" = C8 + C5. 
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with the full theory (at fi = m\y) must be done at the 0{as) level and the renormalization 
group equation has to be solved using the anomalous dimension matrix calculated up to order 
al- Second, the virtual 0{as) corrections for the matrix element (at scale fi ~ rrih) must be 
evaluated and combined with the Bremsstrahlung corrections. The higher order matching has 
been calculated in ref. |T^ and work on the Wilson coefficients is in progress. In this paper we 
will evaluate all the virtual correction beyond those evaluated already in connection with the 
Bremsstrahlung process. We expect them to reduce substantially the strong scale dependence 
of the leading order calculation. 

Among the four-Fermi operators only O2 contributes sizeably and we calculate only its 
virtual corrections to the matrix element for b ^ s'y. The matrix element Oi vanishes because 
of color, and the penguin induced four-Fermi operators O3, ....,06 can be neglected because 





/i = niw 


fx = 10.0 GeV 


^ = 5.0 GeV 


fi = 2.5 GeV 




0.0 


-0.149 


-0.218 


-0.305 


C2 


1.0 


1.059 


1.092 


1.138 




0.0 


0.006 


0.010 


0.014 


C4 


0.0 


-0.016 


-0.023 


-0.031 




0.0 


0.005 


0.007 


0.009 


c. 


0.0 


-0.018 


-0.027 


-0.040 




-0.192 


-0.285 


-0.324 


-0.371 


Cs 


-0.096 


-0.136 


-0.150 


-0.166 




-0.192 


-0.268 


-0.299 


-0.334 


r^ff 


-0.096 


-0.131 


-0.143 


-0.157 



Table 1: Leading logarithmic Wilson coefficients Cj(/i) at the matching scale fi = mw = 80.33 
GeV and at three other scales, /i = 10.0 GeV, /i = 5.0 GeV and fj, = 2.5 GeV. For as{fx) (in the 
MS scheme) we used the one-loop expression with 5 flavors and as{mz) = 0.117. The entries 
correspond to the top quark mass mj{mt^poie) = 170 GeV (equivalent to rrit^poie = 180 GeV). 

their Wilson coefficients Q are much smaller than C2, as illustrated in table |^. However, we 
do take into account the virtual 0{as) corrections to 6 — *• 57 associated with the magnetic 
operators O7 (which has already been calculated in the literature) and Os (which is new). 
Since the corrections to O7 and Os are one- loop diagrams, they are relatively easy to work out. 
In contrast, the corrections to O2, involve two-loop diagrams, since this operator itself only 
contributes at the one-loop level. 

Since the virtual and Bremsstrahlung corrections to the matrix elements are only one (well- 
defined) part of the whole next-to-leading program, we expect that this contribution alone will 
depend on the renormalization scheme used. Even within the modified minimal subtraction 
scheme (MS) used here, we expect that two different "prescriptions" how to treat 75, will lead 
to different answers. Since previous calculations of the Bremsstrahlung diagrams have been 
done in the NDR scheme and also the leading logarithmic Wilson coefficients are available in 

^ This omission will be a source of a slight scheme and scale dependence of the next-to- leading order result. 
^It is consistent to calculate the corrections using the leading logarithmic Wilson coefficients. 
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this scheme, we also use it here. For future checks, however, we also consider in Appendix A 
the corresponding calculation in the HV scheme. 

The remainder of this paper is organized as follows. In section 2 we give the two-loop 
corrections for 6 — > S7 based on the operator O2 together with the counterterm contributions. 
In section 3 the virtual corrections for b ^ s'j based on O7 are reviewed including some of 
the Bremsstrahlung corrections. Then, in section 4 we calculate the one-loop corrections to 
6 — > 57 associated with Og. Section 5 contains the results for the branching ratio for b —>■ sjIq) 
and especially the drastic reduction of the renormalization scale dependence due to the new 
contributions. Appendix A contains the result of the O2 two-loop calculation in the HV scheme 
and, finally, to make the paper self-contained, we include in Appendix B the Bremsstrahlung 
corrections to the operators O2, O7 and Og. 



2 Virtual corrections to O2 in the NDR scheme 

In this section we present the calculation of the matrix element of the operator O2 for b ^ 
up to order in the NDR scheme. The one-loop matrix element vanishes and we must 
consider several two-loop contributions. Since they involve ultraviolet singularities also coun- 
terterm contributions are needed. These are easy to obtain, because the operator renormal- 
ization constants Zij are known with enough accuracy from the order as anomalous dimension 
matrix |Tl[|. Explicitly, we need the contributions of the operators C2SZ2jOj to the matrix 
element for b ^ s'y, where 6Z2j denote the order as contribution of the operator renormaliza- 
tion constants. In the NDR scheme, the non-vanishing counterterms come from the one-loop 
matrix element of C2SZ25O5 and C2SZ2GOQ as well as from the tree level matrix element of the 
operator C2SZ270-/. We also note that there are no contributions to 6 ^ 57 from counterterms 
proportional to evanescent operators multiplying the Wilson coefficient C2. 

2.1 Regularized two-loop contribution of O2 

The dimensionally regularized matrix element M2 of the operator O2 for b sj 

M2 = (S7IO2I6) (2.1) 

can be divided into 4 classes of non-vanishing two-loop diagrams, as shown in Figs. The 
sum of the diagrams in each class (=figure) is gauge invariant. The contributions to the matrix 
element M2 of the individual classes [^-^ are denoted by M2(l), M2(2), M2(3) and M2(4), where 
e.g. M2(l) is 

M2(l) = M2(la) + M2(16) + M2(lc) . (2.2) 



The main steps of the calculation are the following: We first calculate the Fermion loops in 
the individual diagrams, i.e., the 'building blocks' shown in Fig. |^ and in Fig. |^, combining 
together the two diagrams in Fig. ^ As usual, we work in (i = 4 — 2e dimensions; the 
results are presented as integrals over Feynman parameters after integrating over the (shifted) 
loop-momentum. Then we insert these building blocks into the full two-loop diagrams. Using 
the Feynman parametrization again, we calculate the integral over the second loop-momentum. 
As the remaining Feynman parameter integrals contain rather complicated denominators, we 
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(a) 



(b) 



(c) 



Figure 1: Diagrams la, lb and Ic associated with the operator O2. The fermions (6, s and c 
quark) are represented by sohd hues. The wavy (dashed) hne represents the photon (gluon). 






(a) 



(b) 



(c) 



Figure 2: Diagrams 2a, 2b and 2c associated with the operator O2. 



do not evaluate them directly. At this level we also do not expand in the regulator e. The 
heart of our procedure which will be explained more explicitly below, is to represent these 
denominators as complex Mellin-Barnes integrals [0. After inserting this representation and 
interchanging the order of integration, the Feynman parameter integrals are reduced to well- 
known Euler Beta-functions. Finally, the residue theorem allows to write the result of the 
remaining complex integal as the sum over the residues taken at the pole positions of certain 
Beta- and Gamma-functions; this naturally leads to an expansion in the ratio z = {mc/mh)'^, 
which numerically is about z = 0.1. 

We express the diagram in Fig. ^ (denoted by J^) in a way convenient for inserting into the 
two-loop diagrams. As we will use MS subtraction later on, we introduce the renormalization 
scale in the form /i^ exp(7£;)/(47r), where je = 0.577... is the Euler constant. Then, MS 
corresponds to subtracting the poles in e. In the NDR scheme, is given by 

A 



47r2 



r(e) fi^' exp(7£;e) (1 - e) exp(z7re) (^r^if - r%^) L- x 



^The fermion/gluon and the fermion/photon couplings are defined according to the covariant derivative 
D = + igs^A^ + ieQA 
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(a) 



(b) 



Figure 3: Diagrams 3a and 3b associated with the operator O2. We calculate directly their sum 
and denote it by M2(3), see text. 





(a) 



(b) 



Figure 4: Diagrams 4a and 4b associated with the operator O2. We calculate directly their sum 
and denote it by M2(4), see text. 



— x)] 



l-e 



x(l 



+ i5 



(2.3) 



where r is the four-momentum of the (off-shell) gluon, rric is the mass of the charm quark 
propagating in the loop and the term iS is the " e-prescription" . The free index f3 will be 
contracted with the gluon propagator when inserting the building block into the two-loop 
diagrams in Figs. |I| and ^ Note that I 13 is gauge invariant in the sense that r^/^j = 0. 

Next we give the sum of the two diagrams in Fig. using the decomposition in |]18|. We 
get (denoting this building block by Ja/3) 



J, 



^9sQu 

167r2 



E{a, p, r) Az5 + E{a, p, q) A^g - E{(3, r, q) Az23 



(gr) 



-E{a, r, q) -^Ai^^ - E{a, r, q) -^A^26 



9/3 



2 



(2.4) 



(gr) (gr) 

where q is the four-momentum of the photon. The index a in eq. ( p.4|) is understood to be 
contracted with the polarization vector e of the photon, while the index j3 is contracted with 
the gluon propagator in the two-loop diagrams in Figs. ^ and ^. The matrix E in eq. (|2.4|) is 
defined as 

E{a, p, r) = -falf^if - laTp + li3i.ro) - iQafi- (2.5) 

In a four-dimensional context these E quantities can be reduced to expressions involving the 
Levi-Civita tensor, i.e., £'(«,/?, 7) = —iSap-y^i^i^ iy^ the Bjorken-Drell convention). The 
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b s 

XT 



Figure 5: Building block Ip for the diagrams in Figs. ([l|) and (0) with an off-shell gluon. 





Figure 6: Building block Jaii for the diagrams in Figs. (H) and (|). 

dimensionally regularized expressions for the Ai read 

Ais = —4 / dxdy Aiqr)x^ye — Aiqr)xye — 2r'^x'^e + ^ir'^x^e 
Js L 



-r'xe + ?,xC -C 



[{l + e)V{e)eMlEe)l^''C 



2e 



Ai23 = -Ai26 = 8(gr) / dx dy [xye (1 + e)r(e) exp(7£;e)/i^^ C 

•J s 

Az25 = -8(gr) j rfa; [a;(l - a;)e (1 + e)r(e) exp(7£;e)/i^' C 



2e ^-1-el 



2e ^<-l-e^ 



where C and C ^ are given by 

C = ml — 2xy{qr) — x(l — x)r^ — iS 
C-^-' = -exp(i7re) 



^2 ^ ^l?'^) 



1 — X x(l — x) 



i6 



The range of integration in {x,y) is restricted to the simplex S, i.e., < y < [1 
0<x<l. 



(2.6) 



Aig = 4 / dxdy Aiqr)xy'^6 — Aiqr)xye — 2r'^x'^ye + 2r'^x'^e + r^xye 
Js L 

-2r2xe + 3|/C - d [(1 + e)r(e) exp(7£;e)/x2^ C"^"^] (2.7) 



(2.8) 
(2.9) 



(2.10) 
x) and 



Due to Ward identities, not all the Ai are independent. The identities given in ref. [18 
the context of the full theory simplify in our case as follows: 



m 



g"J,^ = ; r^Jap = Q . (2.11) 

They allow to express Ai^ and Aig in terms of the other Ai which have a more compact form. 
These relations read 

Air, = Az23 , 

2 

Ale = 7^A«25 + Az26 • (2.12) 
(gr) 

Of course, eq. ( p.l2| ) can be checked explicitly for all values of e, using partial integration and 
certain symmetry properties of the integrand. 

We are now ready to evaluate the two-loop diagrams. As both and Jap are transverse 
with respect to the gluon, the gauge of the gluon propagator is irrelevant. Also, due to the 
absence of extra singularities in the limit of vanishing strange quark mass, we set = from 
the very beginning (the question of charm quark mass "singularities" will be discussed later). 

As an example, we present the calculation of the two-loop diagram in Fig. Ic in some detail. 
Using 1/3 in eq. ( |2.3| ), the matrix element reads 

{2nY 



M(lc) = ^eg,(?2^pr(e)exp(27se)/i^^(l-e)exp(me)(47r)-^y' 



Kl-^)]^- . (2.13) 

[r^ — m^/ (x(l — x)) + i5Y 

In eq. ( p.l3| ), u{j)') and u{p) are the Dirac spinors for the s and the h quarks, respectively, 
and Cp = 4/3. In the next step, the four propagator factors in the denominator are Feynman 
parametrized as 

1 r(3 + e) r dudvdwdy y"^^ 5{\ — u — V — w — y) (o^A\ 

DiD2D^Dl " r(e) J [Diu + D^v + D^w + D^yf^' ' ^' ' 

where Di = {p' — r)^, D2 = {p — r)^, D3 = and = — m^/(x(l — x)). Then the integral 
over the loop momentum r is performed. Making use of the 5 function in eq. ( [^.14D , the integral 
over w is easy. The remaining variables v and y are transformed into new variables m', v' 
and y\ all of them varying in the interval [0,1]. The substitution reads 

u^v\\-v!) , v^v'v! , y^{\-v')y' . (2.15) 

Taking into account the corresponding jacobian and omitting the primes(') of the integration 
variables this leads to 

M(lc) = — ^ eQd5f^Cirr(2e) exp(27£;e)/i^' exp(2z7re) (1 - e) / dxdudvdyx 
647r* J 

[xil - x)Y'' y'-^ (1 - vY V 



u{p' 



C 1 1 



u{p) , (2.16) 



8 



where Pi, P2 and P3 are matrices in Dirac space depending on the Feynman parameters x, u, 
V, y in a polynomial way. C is given by 



C = mlv{l — v)u — 



X 1 



x] 



-(1 — v)y + i6 



(2.17) 



In what follows, the ultraviolet e regulator remains a fixed, small positive number. 

The central point of our procedure is to use now the Mellin-Barnes representation of the 
"propagator" l/(fc2 - M^)^ [|l|, |g, |l|, |2| which is given by 



(fc2-M2)^ (F)^ r(A) 



f ds{-MykyT{-s)T{X + s) 



(2.18) 



where A > and 7 denotes the integration path which is parallel to the imaginary axis (in 
the complex s-plane) hitting the real axis somewhere between —A and 0. In this formula, 
the "momentum squared" k"^ is understood to have a small positive imaginary part. In ref. 



19|, |21| exact solutions to Feynman integrals containing massive propagators are obtained by 



representing their denominators according to the formula ( p.l8| ) with subsequent calculation of 
the corresponding massless integrals. 

In our approach, we use formula ( p.l8| ) in order to simplify the remaining Feynman pa- 
rameter integrals in eq. (|2.16|) . We represent the factors l/C*^'^ and l/C^+^e gq_ (|2.16|) as 
Mellin-Barnes integrals using the identifications 



1 — v)u 



x{l 



x] 



v)y 



(2.19) 



By interchanging the order of integration, we first carry out the integrals over the Feynman 
parameters for any given fixed value of s on the integration path 7. These integrals are basically 
the same as for the massless case mc = (in eqs. ( p.l6| ) and ( |2.17| )) up to a factor (in the 
integrand) of 



y 



uvx{l — x) 



exphiTsj — ^ 
\ mt . 



(2.20) 



Note that the polynomials Pi, P2 and P3 have such a form that the Feynman parameter integrals 
exist in the limit 0. If the integration path 7 is chosen close enough to the imaginary 

axis, the factor in eq. ( |2.20D does not change the convergence properties of the integrals, i.e., 
the Feynman parameter integrals exist for all values of s lying on 7. It is easy to see that the 
only integrals involved are of the type 



dww^ 



or 



dwwP{l - wy = p{p+l,q + 1) 



r(p + i)r(g + i) 

T{p + q + 2) 



(2.21) 



For the s integration we use the residue theorem after closing the integration path in the right 
s-halfplane. One has to show that the integral over the half-circle vanishes if its radius goes 
to 00. As we explicitly checked, this is indeed the case for (m^/m^) < 1/4, which is certainly 
satisfied in our application. The poles which lie inside the integration contour are located at 



0,1,2,3,4, 



l-e,2-e,3-e,4 
1 - 2e,2 - 2e,3 



5 ^ ^5 

2e,4 



2e, 



(2.22) 
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The other two-loop diagrams are evaluated similarly. The non-trivial Feynman integrals 
can always be reduced to those given in eq. ( p.21j ) after some suitable substitutions. The only 
change is that there are poles in addition to those given in eq. ( |2.22| ) in those diagrams where 
the gluon hits the 6-quark line; they are located at 

s = 1/2 - 2e, s = 3/2 - 2e, s = 5/2 - 2e, s = 7/2 - 2e, (2.23) 

The sum over the residues naturally leads to an expansion in z = {ml /ml) through the 
factor [ml/mlY in eq. (|2.2CI|). This expansion, however, is not a Taylor series because it also 



involves logarithms of z, which are generated by the expansion in e. A generic diagram which 
we denote by G has then the form 



2 

G = co + Ecn™^"log"^^ , , (2.24) 



mt 



where the coefficients Cq and Cnm are independent of z. The power n in eq. ( p^.24D is in general 
a natural multiple of 1/2 and m is a natural number including 0. In the explicit calculation, 
the lowest n turns out to be n = 1. This implies the important fact that the limit mc 
exists; thus, there cannot be large logarithms (from a small up-quark mass) in these diagrams. 

From the structure of the poles one can see that the power m of the logarithm is bounded 
by 4 independent of the value of n. To illustrate this, take n = 100 as an example. There are 
3 poles located near n=100, viz., at s = 100, s = 100 — e,s = 100 — 2e, respectively (see eq. 
(|2.22| )). Taking the residue at one of them, yields a term proportional to 1/e^ coming from the 
remaining two poles. In addition there can be an explicit 1/e^ term from the integration of the 
two loop momenta. Therefore the most singular term can be l/e"^. Multiplying this with z^ in 
eq. (|2.20|) leads to log™" z where m can be 4 at most. 

We have retained all terms up to = 3. Comparing the n = 3 numerical result with the 
one obtained by truncating at n = 2 leads to a difference of about 1% only. 

We have made further checks of our procedure. For example we have calculated diagram lb 
directly. Expanding the result, we reproduce the expressions obtained by applying the Mellin- 
Barnes integral at the Feynman parameter level as described above. A similar exercise for the 
imaginary part of diagram Ic shows that the exact and the expanded result (up to z^ terms) in 
these examples agree at the 1% level. In addition, we checked that the imaginary part of the 
sum of all diagrams coincides numerically with the results of Soares [^, |^ (note, however that 
in the physical region only the diagrams in Fig. and Fig. ^ contribute to the imaginary part). 
In ref. Soares applied dispersion techniques to calculate the real part. However, using the 
imaginary part in the physical region, only the real part of the diagrams in Fig. |1| and Fig. ^ 
is obtained. We have checked that our numbers for these two sets of diagrams indeed coincide 
with the results of Soares. However, the contribution of the diagrams in Fig. ^ and Fig. ^ is 
missing in p| because the additional unphysical cuts were not taken into account. We also note 
that a separate consideration of the subtraction terms would be required to obtain the correct 
/i dependence. 

We mention that the Dirac algebra has been done with the algebraic program REDUCE ^ 



24|. The Feynman parameter integrals and the determination of the residues have been done 



with the symbolic program MAPLE 

^Some checks have been done with TRACER 
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We now give the results for the diagrams shown in Figs. 1 to 4. As mentioned already, the 
individual diagrams in each figure are not gauge invariant but only their sum. Note, that the 
leading ultraviolet singularity in the individual two-loop diagrams is in general of order 
In the gauge invariant sums M2{i) the cancel and we are left with 1/e-poles only. The 
results read (using z = {rric/mbY and L = logz): 



-4e 



[37 - (540 + 216L)z 



1 rrib 
36e V ; ' 216 

+ (2167r^ - 540 + 216L - 216L^)z^ 



-{lUn^ + 136 + 240L - 1UL'^)z^ 



+ 



ITT 

18 

a. 



1 - 18z + (18 - 36L)z2 ^ ^2AL - 20)r 



X—CFQd{si\07\h)tree 
TT 



(2.25) 



M2(2) 



-At 



— + 



1^ 



13 + (367r2 - 108); 



36e V A* y 216 
-14471^^=^/^ + (648 - 648L + 108L'^)z'^ 
+ (1207r^ - 314 + 12L + 288L'^)z^] } 



y<-CFQa{si\0,\h) 

TT 



tree 



(2.26) 



-At 



+(72 - 127r^ - 96C(3) + (96 - 247r^)L + 12L^ + 8L^)z 
+(60 + 247r^ - 96C(3) + (24 - 247r^)L - 2AL^ + 8L^)z^ 
-(68 - A8L)z 

in 
+T2 



-3 + (24 - 27r^ + 6L + QL^)z 
+(6 - 27r2 - 12L + 6L2)^2 ^ ^^z^-^ } 



><—CFQu{si\07\h)tree 
TT 



(2.27) 



+(-24 + 271^ + 24C(3) - (12 - Qt:'^)L + 2L^)z 
+ (-12 - 47r^ - 48C(3) + 12L - 6L^ + 2L^)z^ 
+ (-6 + Gtt^ - 24L + l8L'^)z^] } 

Oi 

X—CFQu{s-i\07\h)tree 
TT 



(2.28) 
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In these expressions, the symbol ( denotes the Riemann Zeta function, with ({3) ~ 1.2021; 
Qu = 2/3 and Qd = —1/3 are the charge factors for up- and down- type quarks, respectively. 
The matrix element {s'y\Of\b)tree is the 0(a°) tree level matrix element of the operator O7; its 
explicit form is 

{sj\Or\b)tree = m,^u{p')^(iRu{p) . (2.29) 
In formula ( |2.29 ) mt, should be identified with the running mass mb{fJ>) in principle (see eq. 



(|1.2|) ). However, as the corrections to O2 are explicitly proportional to a^, can be identified 
with the pole mass as well (apart from 0{al) corrections which we systematically neglect). 

2.2 Count erterms 

The operators mix under renormalization and thus the counterterm contributions must be 
taken into account. As we are interested in this section in contributions to 6 — >■ S7 which are 
proportional to C2, we have to include, in addition to the two- loop matrix elements of C2O2, 
also the one-loop matrix elements of the four Fermi operators C2SZ2jOj {j = 1, ...,6) and the 
tree level contribution of the magnetic operator C2SZ27O7. In the NDR scheme the only non- 
vanishing contributions to b —>■ sj come from j = 5, 6, 7. (For j = 5, 6 the contribution comes 
from the diagram in which the internal 6-quark emits the photon). The operator renormalization 
constants Zij are listed in the literature [|lT| in the context of the leading order anomalous 
dimension matrix. The entries needed in our calculation are 

cry ^ c-ry ^ 

0^25 — —JZ '-'F ' ^^26 — 77 L'F , 

487re lovre 
fy 8 

5^27 = —^{QQu--Qd)Cp . (2.30) 

lovre 9 

Defining 

M2, = {s^\6Z2Pj\b) , (2.31) 
we find the following contributions to the matrix elements 

M25 = -^QdCF-( — ] {Sl\0r\b),ree 

M26 = ^QdC^-f — 1 '\s^\Or\b),ree 
lOTT e \ jJ, J 

M,r = ^ - 9j^) 1 {n|0,|(,>.„. (2,32) 

We note that there is no one-loop contribution to the matrix element for 6 ^ 57 from the 
counterterm proportional to C2 \ OI2, where the evanescent operator (see e.g. the last ref. 
in 0]) reads 

0r2 = l02{i,^7[,7ulp])-02 . (2.33) 
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2.3 Renormalized contribution proportional to C2 



Adding the two- loop diagrams from section 2.1 (eqs. (|2.25|) , (|2.26|) , (|2.27|) and ( p.28|) ) and the 
counterterms from section 2.2 (eq. ( |2.32D ) we find the renormahzed contribution M2 

M2 = M2(1) + M2(2) + M2(3) + M2(4) + M25 + M26 + M27 . (2.34) 

Of course, the ultraviolet singularities cancel in M2. Inserting Cp = 4/3, Qu = 2/3 and 
= — 1/3, we get the main result of this paper, which in the NDR scheme reads: 



with 



M2 = {S^\07\b)tre 



?2 log h r2 

/i , 



416 



(2.35) 
(2.36) 



243 
+ 

+ 



2 3/2 



{-833 + 1447r 

1728 - ISOn^ - 1296C(3) + (1296 - 3247r^)L + 108^^ + 36L' 
648 + 7271^ + (432 - 216n^)L + 36L' 
-54 - 8471^ + 1092L - 756L^ 



(2.37) 



1^1-5 
81 I 



45 - Sn^ + 9L + 9L^ 



z + 



-Sn^ + 9^2 



[28-12L] 2^} (2.38) 



Here, 3ftr2 and denote the real and the imaginary part of r2, respectively. The quantity z 
is defined a.s z = {m'^/ml) and L = log(2;). In Fig. ^ we show the real and the imaginary part 
of r2. For z > 1/4 the imaginary part must vanish exactly; indeed we see from Fig. |^ that the 




0.05 0.1 0.15 0.2 0.25 
z = (mymj2 



Figure 7: Real and imginary part of r2 in the NDR scheme (from eqs. ( p.37|) and (|2.38|) ) 



imaginary part based on the expansion retaining terms up to z^ indeed vanishes at 2; = 1/4 to 
high accuracy. 
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3 0{as) corrections to O7 



The virtual corrections associated with the operator O7 as shown in Fig. §3 (together with the 
selfenergy diagrams and the counterterms) have been taken into account in the work of Ali and 
Greub, see e.g. & H, O], where 7^ was retained. Since we neglect tjIs in this work, we are 



7 



s 



0, 

(a) 




Figure 8: Virtual corrections to O7 



interested only in the limit 0. Because of the mass singularities in the virtual corrections 

(which will be cancelled when also taking into account Bremsstrahlung corrections), we only 
keep TUs as a regulator. 

Including the lowest order contribution, the result then becomes (using p = {nis/mhY) in 
the NDR scheme 




{s-i\0^\h)^irt 



{S'y\0j\h)tree 



2 4 
logV logp - logp 



8 + 4 log 



p ■ 



(3.1) 
(3.2) 



JJ) contains all the counterterm contributions. The 1 / e/ij-poles in this equation 
are therefore of infrared origin as indicated by the notation. The last term in the curly bracket 
in eq. ( |3.2| ) represents a p dependence of ultraviolet origin. The additional p dependence, which 
is generated when expanding in em is cancelled at the level of the decay width together with 
the l/e/i{-poles when adding the Bremsstrahlung correction due to the square of the diagrams 
associated with the operator O7. As all intermediate formulae are given in the literature, we 
only give the final result for the corrections (virtual+Bremsstrahlung) to the decay width. 
Denoting this contribution by we get in the limit = 

16 47r^ , mu 
— +41og- 



r 



77 



1 + 



37r V 3 



where the lowest order contribution F^^ reads 



F° 



m2(p)mf 



32^4 I ^ t 7 I 



Op 



(3.3) 



(3.4) 



For later convenience, we can formally rewrite {s^\Oi\b)mrt in eq. ( p.l|) in such a way 
that its square reproduces the result in eq. (|3.3| )- This modified matrix element, denoted by 
{s^\0^\b)rad then reads 



{s-i\0^\h)rad = {s-f\07\b)tree 1 + 



as 
47r 



log 



(3.5) 



14 



with 



r, = ^(4-vr^) 



(3.6) 



4 Virtual corrections to Og 

Finally, we consider the contributions to b s'j generated by the operator Os, i. e., 

Ms = {s-f\Os\b) . (4.1) 
The corresponding Feynman diagrams are shown in Fig. 0. The sum of the three diagrams in 







8 







8 



^ • X ,X X, X • 5^ 



(a) 



(b) 



Figure 9: Contributions of Og to 6 57. The cross (x) denotes the possible place where the 
photon is emitted. Figs, (a) and (b) are separately gauge invariant. 



Fig. pk yields 



Ms (9a) 



2 TT 

while the diagrams in Fig. ^3 gives 



2 + 2 log(mf,//i) — in 

e 



{si\0,\b) 



tree ; 



12 TT 



6 

21 + 27r2 + 121og(mfe//i) 



{S'j\07\b)tree 



Thus the sum of all the six diagrams is 



Ms (9) 



QdCp as 



12 

— - 33 + 27r2 + 24 log(mfe//i) - Qin 



12 TT 

There is also a contribution from a counterterm; it reads 



{si\Or\b) 



tree 



The renormalization constant 



Msr = 6Zs7 {Sj\07\b)tr, 



SZgj — — CpQd - 
TT e 



(4.2) 



(4.3) 



(4.4) 



(4.5) 



(4.6) 



has been calculated in the literature . The sum of all contributions leads to the renormalized 
result Ms: 



Ms = (s7|07|&)i.ee^ (^fslogy + rs 



with 



32 



4 / 9 

(-33 + 27r2-6m 

27 



(4.7) 
(4.8) 
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5 Results and conclusions 



We have calculated the virtual corrections to 6 ^ 57 coming from the operators O2, O7 and 
Og. The contributions from the other four- Fermi operators in eq. ( |1.2| ) which are given by the 
analogous diagrams as shown in Figs. (|l|) - (^) were neglected, because they either vanish (Oi) 
or have Wilson coefficients which are about fifty times smaller than that of O2 while their matrix 
elements can be enhanced at most by color factors. However, we did include the non-vanishing 
diagrams of O5 and Og where the gluon connects the external quark lines and the photon is 
radiated from the charm quark because these corrections are automatically considered when 
C'^^^ (defined in eq. (|1.3| )) is used instead of C-j. As discussed in section three, some of the 
Bremsstrahlung corrections to the operator 0^ have been transferred into the matrix element 
for 6 — *• S7 in order to present results which are free from infrared and collinear singularities. 

The sum of the various contributions derived in the previous sections yields the amplitude 
A{b 57) for h — i> S7. The result can be presented in a convenient way, following the treatment 
of Buras et al. ||^, where the general structure of the next-to- leading order result is discussed 
in detail. We write 

A{h Sj) = D {S^\0r{fl)\b),ree (5.1) 



with D 



/^ = C^r^^(/^) + ^E(^P''(^Kaog^ + Cf)^^^r.j , (5.2) 

and where the quantities ii and are given for z = 2, 7, 8 in sections 2, 3 and 4, respectively. 
For the full next-to-leading logarithmic result one would need the first term on the rhs of eq. 
( |5.2| ) jCj^^ln), at next-to-leading logarithmic precision. In contrast, it is consistent to use the 
leading logarithmic values for the other Wilson coefficients in eq. ( |5.2| ). As the next-to-leading 
coefficient C^^^ is not known yet, we replace it by its leading logarithmic value cf^^^^ in the 
numerical investigations. The notation {s^\0-r[ix)\b)tree in eq. ( |5.1| ) indicates that the explicit 
factor rrtb in the operator Oj is the running mass taken at the scale /x. 

As the relevant scale for a h quark decay is expected to be /i ~ m^, we expand the matrix 
elements of the operators around /i = m^, up to order 0(0;^). Thus we arrive at 

A{h^s^) = -I^ D {s^\Orimt)\b)tree (5.3) 

D = C,^//(/.) + ^ E [cr''"i^^hP'" log ^ + Cr'-^r^l , (5.4) 
with the quantities 'ylj'^^^^ 

^p-ff = i^ + 85„ (5.5) 

being just the entries of the (effective) leading order anomalous dimension matrix p. As also 
pointed out in this reference, the explicit logarithms of the form ^^(m;,) log(m;,//i) in eq. ( |5.4|) 
should be cancelled by the /i-dependence of Cf^^^^jj). This is the crucial point why the scale 
dependence is reduced significantly as we will see later. 

^As we neglect the virtual correction of Oa-Oe, there is a small left-over /i dependence, of course. 



with D 
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From A{b 57) in eq. ( p.3|) we obtain the decay width F"*''* to be 



■^virt 



F\D\ 



(5.6) 



where we discard term of 0{al) in \D\ . The factor F in eq. (|5.6| ) is 



^ ^b,pole 



1 - - 

3 



TT 



(5.7) 



To obtain the inclusive rate for B Xg'y consistently at the next-to-leading order level, we 
have to take into account all the Bremsstrahlung contributions. They have been calculated by 
Ali and Greub for the operators O2 and O7 some time ago 0] , while the complete set has been 
worked out only recently [l^, 0. Here, we neglect the small contribution of the operators 
O3 - Og as we did for the virtual corrections, i.e., we only consider O2, Oj and Og. The 
corresponding Bremsstrahlung formulae are collected in Appendix B. 

In order to arrive at the branching ratio BR(6 — 37(5')), we divide, as usual, the decay width 
F(5 — i> Xs'~f) = F"*** + r^^'^'^^ by the theoretical expression for the semileptonic decay width F^^ 
and multiply this ratio with the measured semileptonic branching ratio BR^^ = (10.4 ± 0.4)% 
13, i.e., 

BR{b^s-f{g)) = ^BRsi . (5.8) 

J- si 

The semileptonic decay width is 



2 4 



<^ 3 - 



X 

I 
I 

cn ^ 
m 
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Figure 10: Branching ratio for 6 ^ S7 based on the leading logarithmic formula in eq. ( |5.12| ). 
The upper (lower) solid curve is for fi = nib/ 2 {ft = ^rrib). The dotted curves show the CLEO 
bounds [0. 



= '^"^''"^'^ V - 



f{mc/mb) 



(5.9) 
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Figure 11: Branching ratio for b s'y{g) neglecting the virtual corrections of O2 and Og 
calculated in the present paper. The upper (lower) solid curve is for /i = mb/2 (/i = 2mfe). The 
dotted curves show the CLEO bounds 0. 

where the phase space function g{u) is defined as 

g{u) = l-8u^ + 8u^ -u^ -2AuHogu , (5.10) 
and an approximate analytic form for the radiative correction function f{u) has been found in 



281 to be 



/M=(-^-|) (1--)^ + ^ . (5.11) 

In Figs. 10-12 we compare the available leading- log results and our new results for the 
inclusive branching ratio for B Xgj as a function of the top quark mass. A rather crucial 
parameter is the ratio rric/mh; it enters both, b ^ s'j through the virtual corrections of O2 and 
the semileptonic decay width through phase space. It can be written as mc/mb = 1 — (m^ — 
^i^/^b- While the mass difference — mc is determined quite precisely through the l/mq 



expansion ||29| or from the semileptonic b ^ c spectrum; we use rrib — rric = 3.40 GeV ||30[| , 
the 6— quark mass is not precisely known. Using for the b quark pole mass mb^poie = 4.8 ± 0.15 
GeV one arrives at rric/mb = 0.29 ± 0.02. In the plots the central values for rrib^poie and rric/mb 
have been used. Moreover, we put Vcb = Vts and Vtb = 1 and also use the central value for 
the measured semileptonic branching ratio BR^^ = 10.4% in eq. ( ^.8|) . In all three plots the 
horizontal dotted curves show the CLEO la-limits for the branching ratio BR(i? X^'y) ^ . 

In Fig. |1^ we show the leading logarithmic result for the branching ratio of 6 57, based 
on the formula 

BR(6 S7)'^"'^^"3 = .^"'^^ , BR,^ . (5.12) 

TTg{mJmb) 



Similarly, Fig. 11 exhibits the results also taking into account the Bremsstrahlung corrections 
and the virtual corrections to O7 without including the virtual corrections of O2 and Og. We 
can reproduce this result by putting £2 = ''"2 = = ''^8 = in our formulae. As noticed in the 
literature P, |15|, the /i dependence in this case is even larger than in the leading logarithmic 
result shown in Fig. |1^. Also, the experimentally allowed region is shown. 
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Figure 12: Branching ratio for b S'~f{g) based on the complete formulae presented in this 
section. The upper (lower) solid curve is for fj, = mb/2 {[i = 2mb). The dotted curves show the 
CLEO bounds 0. 



The theoretical results shown in Figs. ^ and |TT] allowed for a reasonable prediction of the 
branching ratio within a large error which was essentially determined by the /i dependence. As 
we see, they do agree well with experiment but within large uncertainties. 

In Fig. |l^, finally, we give the branching ratio based on formula (|5.8|) which includes all 
virtual corrections calculated in the present paper. Because all the logarithms of the form 
as{mb)log{mh/ fi) cancel as discussed above, the fi dependence is significantly reduced in our 



improved calculation (Fig. |1^); the bands of scale uncertainty are rather narrow. To illustrate 
the remaining renormalization scale dependence, we present two 'scenarios' which differ by 
higher order contributions. First, we take the explicit factors in eqs. ( |5.4| ) and ( p.9[ ) at 
fi = rrif, as indicated in these formulae; varying /i between {mb/2) and {2mb) in formula (|5.4|) 
leads to the dash-dotted curves in Fig. |T^. Second, we evaluate the explicit in eqs. ( |5.4| ) 
and i f) .9]) at the (variable) scale fi. Varying again the scale fj, between {mb/2) and {2mb) yields 
the solid lines in Fig. |l^. In both scenarios the upper (lower) curve corresponds to fi = mb/2 
{n = 2mb)- We mention that the /i-band is larger in the second scenario and it is therefore 
safer to use this band to obtain a feeling for the remaining scale uncertainties. 

While this result shows that the theoretical accuracy can be strongly improved by the 
next-to-leading calculations, it would be premature to extract a prediction for the branching 
ratio from Fig. |T2| (with obviously a small error) and claim, for instance, a discrepancy with 
experiment. This only will become possible (with high theoretical precision) if also C^^^ is 
known to next-to-leading logarithmic precision. This additional effect will , essentially, shift 



the narrow bands in Fig. |T2|, without broadening them substantially. The drastic reduction of 
the theoretical uncertainties shows that significant experimental improvements are necessary 
to extract the important information available in these decays. 
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Appendix 



A O2 contribution in the 't Hooft-Veltman scheme 

In this appendix we present the results for the matrix elements of the process 6 — >^ 57 based on 
the operator O2. In addition to the two-loop diagrams shown in Figs. |I| ^ 1^, we also give all 
the counterterm contributions which multiply the Wilson coefficient C2. 



A.l The 't Hooft-Veltman Scheme 



In the 't Hooft-Veltman scheme |T^, the d dimensions are split into 4 and d — A; the corre- 



sponding structures are distinguished by no superscript, by a tilde and by a hat, respectively. 
There are Lorentz indices in d, 4 and d — 4 dimensions and the corresponding metric tensors 
d/iu, Qiiv and g^^y. While all the 7-matrices are taken in d dimensions, their indices are split in 
4 and d — A components, according to the rules 

-g^.r" = 4 , g^,g'^'' = d-4: , 

-g.ug"' = . (A.l) 

The 7-matrices in 4 dimensions (7^) and {d — 4) dimensions (7'') are defined by g^^^^v and 
g^^'^ui respectively. Assuming the usual anticommutation relations of the ci-dimensional Dirac 
matrices in terms of the (i-dimensional metric tensor g^y, one gets the following rules for 7^^ 
and 7^^. 

{7^r} = 2^'^^ , {^,7^ = 2^^^^ , {7^r} = o . (a.2) 

The commutation relations with 75 are postulated to be 

{7^75} = , [r,75]=0 , (A.3) 

which is equivalent to defining 75 by the product i7o7i7273- This is the only way known to 
treat 75 without running into algebraic inconsistencies |^T|. Finally, we mention that the chiral 
vertices in d-dimensions can be defined in different ways, all having the same formal limit when 
(i — 4. For left- and right- handed currents we follow the common practice and use 

YL = RYL and ^^'R = L-f''R . (A.4) 

There are several possibilities to define the operators O7 and Og in d dimensions (with identical 



4 dimensional limit); for example the term a'^'^ in eq. ( |1.2|) can be defined to be 



^"^^ = ^[7^71 , (A.5) 
where 7^^ and 'j" are the "(i-dimensional" matrices, or, alternatively, 

^^^^ = ^[7^71 • (A.6) 

A difference in the definition will result in a difference of the finite terms in the one-loop matrix 
elements of these operators. 

As the calculation of the one-loop matrix elements of O7 and Og is relatively easy, once an 
exact definition of the operators has been specified, we give in this appendix only the result for 
the two- loop contribution of the operator O2. 
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A. 2 One-loop building blocks 

The principal steps of the calculation in the HV scheme of the O2 contribution to the matrix 
element for b ^ s'-f are the same as in the NDR scheme. Therefore, we again start with the 
one-loop building blocks. For the building block I/3 in Fig. § we get in the HV scheme 



9s 



in- 



r(e) /i^" exp(7£;e) exp(i7re) / dx\x{l — 

Jo 



x(l 



+ i6 



x{l - x) {r^r^-f^ - r^-fp) + -^r' (7/377 + 
while the building block J^/s in Fig. ^ reads 



(A.7) 



327r2 (1 + e) 



E{a, (3, r) At^ + E{a, /?, q) Aiq - E{(3, r, q) -^A^ss 



(qr) 



-E{a, r, q) -^Aiss - E{a, r, q) T^Aise 



-nir 



(qr) 
(qr) 



(qr) 

7/3^^0 



7a 



' 2 



7c 



+7a7/3(gr)} TL-A 



(A.8) 



The quantities E and Ai are given in eq. ( |2.5|) and (2.6)- (|2.9| ), respectively. As the Ai are the 
same as in the NDR scheme, the Ward identities in eq. ( |2.12| ) still hold. The additional term 
A in eq. (|A.8| ) reads 

(A.9) 



A = / dxdy(qr)C~^-'er(e)fi^'e"'^' 
Js 



where C~^~^ and the integration range S are defined in eq. ( [2.1(J ). 

We note that the terms linear in in eqs. ( [A.7| ) and ( |A.8| ) are absent in schemes where 
75 has formally the same anticommuting properties as in the 4-dimensional Dirac algebra. 
Especially, they are absent in the NDR scheme. However, working consistently in the HV 
scheme, we have to retain these extra terms. 

In order to make the comparison with the NDR calculation easier, we discard in a first step 
these extra rric terms and denote again the corresponding results by M2{i) where i = 1,2,3,4 
numbers the figure in which the corresponding diagrams are shown. In sections A. 3, A. 4 and 
A. 5 the terms linear in rric are left out; they are denoted by M|(2) and collected in section A. 6. 
Those counterterm contributions involving evanescent operator which generate only terms linear 
in rric are also discussed there. 



A. 3 Regularized two-loop contribution of O2 

We now give the results for the two-loop diagrams in Figs. § discarding the linear terms in rric 
in the one- loop building blocks in eqs. ( [A.7| ) and ( |A.8| ). The results read, using z = {rric/mf,)^ 
and L = log z, 



Moil] 




-4e 



+ 



216 



[40 - (540 + 216L)z 



22 



+(21677^ - 540 + 216L - 216L'^)z'^ 
-(1447r^ + 136 + 240L - UAL^)z^ 



+ 



18 



1 - 18^ + (18 - 36L)^2 + (24L - 20)^' 



^—CFQd{si\07\h)tree 
TT 



(A.IO) 



M2(2) 



1 

18e 



-4e 



-17 + (367r^ - 108)^ 



. /i ; 216 

-1447r2//2 + (648 - 648L + l{)^L^)z^ 
+ (1207r2 - 314 + 12L + 288L2)^^] } 

X—CFQd{si\07\h)tree 
TT 



(A.11) 



+ (72 - 127r=^ - 96C(3) + (96 - 2471^)^ + 12L^ + 8L^)z 
+ (60 + 247r2 - 96C(3) + (24 - 247r^)L - 24L2 + 8L^)z'^ 
-(68 - 48L)z3l 



ITT 

+T2 



3 + (24 - 27r^ + 6L + 6L^)z 
+ (6 - - 12L + 6L2)^2 + 12z^] } 



cx 

X—CFQu{s'y\07\b)tree 
TT 



(A.12) 



Ms (4) 




-4e 



+(-24 + 27r^ + 24C(3) - (12 - 67r^)L + 2L^)z 

+(-12 - 47r^ - 48C(3) + 12L - 6L^ + 2L^)z^ 

+(-6 + dn^ - 24L + 18L2)r^] } 
cy 

TT 



X 



I tree 



(A.13) 



A comparison with the corresponding NDR scheme expressions in section 2.1 shows that the 
imaginary parts are identical for each set M2{i). This property has to be satisfied because these 

imaginary parts can be derived by cutting rule techniques where no rcgularization is necessary. 
As expected, the only difi^erence between the two schemes is in the independent terms of 
the real part. 
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A. 4 Count erterms 



There are several counterterm contributions (multiplying C2) which all have to be taken into 
account. As noted, in this section we only give the counterterms whose effects do not lead to 
terms linear in rric- There is just one such counterterm coming from 5Z27O7. The operator 
renormalization constant Z27 reads 



6Z27 = — {^Qu + —Qd) Cp ■ 
Tie 8 36 

Defining 

the counterterm contribution is given by 

as fSQuCp QdCF\ 1 , 1^ |,\ 



(A.14) 



(A.15) 



(A.16) 



A. 5 Renormalized contribution proportional to C2 

Adding the two-loop diagrams (eqs. (lAlOl) , ([AJID , ( |A1^ ) and (lAlsD ) and the counterterm 
(eq. ( [A.16| )), we arrive at the renormalized contribution which we denote by M2. 



M2 = M2(1) + M2(2) + M2(3) + M2(4) + M27 • 
Inserting Cp = 4/3, Qu = 2/3 and Qd = —1/3 we get in the HV scheme 



(A.17) 



with 



M2 = {s-f\07\b)tre 



as 
47r 



?2 log h r2 

/i , 



416 



(A.18) 



(A.19) 



243 
+ 



{-860 + lUn^z^/'^ 

1728 - 1807r2 - 1296C(3) + (1296 - 3247i^)L + 1081"^ + 36L^ 
648 + 7271^ + (432 - 216n^)L + 36L' 



+ 



-54 - 847r^ + 1092L - 756L^ 



} 



(A.20) 



16n ( 



45 - 37r2 + 9L + QL^] z + + ^L^\ + [28 - 12L] z^ \k.2l) 



Again, 3fJr2 and 3r2 denote the real and the imaginary part of r2, respectively. Comparing with 
the final renormalized expression M2 in the NDR scheme, given in section 2.3, we conclude 
again that the scheme dependence only affects the independent term in the real part of r2. 
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A. 6 Terms linear in rric 

Finally, we collect the linear rric- terms of the two-loop digrams which stems from in the building 
blocks in eqs. ( |A.7D and ( |A.8| ). We obtain 



M|(l) = 

2V ; 24 



2 , rrih 

- +17-81og— + 47ri 

e /i 



a 



-CFQ,{s^\0:}r\b)tree 



TC 



M|(2) = 

2V ; 24 



e 



I log — 



-C^gd(s7|07TI&)tree 



TT 



M|(3) + M2^(4) = - 



4 



2 , 

- + 3 - 8 log — 

e 



TT 



Ci.Q.(s7|0r"|&)tr-ee 



where the new-induced magnetic type operator O"^*" in eq. ( |A.24| ) reads 

e 



167r2 



(A.22) 

(A.23) 
(A.24) 

(A.25) 



The operators O^f^ or O^fl", which only enter in intermediate steps, contain an additional left- 
or right- handed projection operator after the a^jy term in eq. ( |A.25| ). 

In addition to these two-loop diagrams, there are two counterterms proportional to evanes- 
cent operators which lead to terms linear in rric. The first (denoted by ^/^Ofpcrmi) comes from 
one-loop gluon corrections to the four-Fermi operator O2 and is given explicitly in ref. ||32 



As we use this operator only as an insertion into the matrix element for b 57, we adapt the 
general color structure given in |^ for this special case. In our notation this counterterm takes 
the form 

^CpEhv , (A.26) 



^ 4: Fermi 



57r e 



with 



-12e (r_ ® r+ + r+ ® r_) - [7r7p7M75 ® 7"7''7''75 + 27^ 27^75 ® 7''75 ] 

- [7^7pr_7 V ® r_ - 7^7^r_ ® r_7^7^ - r_7p7^ ® 7^7''r_ + r_ ® 7^7^1.7 v ] , 

(A.27) 

where r± = 7^ (1 ± 75)/2. Its contribution to the amplitude for 6 ^ 57 is 



rrir 



a. 



Mf [four-Fermi] = - - 21og^ -CpQ^ {si\0!}nb)tree 



7r 



(A.28) 



The second counterterm (denoted by 1/e corresponds to the 1/e pole term in rric- 

term in the one- loop building block Ip in eq. ( [A.7|) . It can be written as 



1 



Qev 

penguin 



1 g \^ 
e 167r2 ~ 



(A.29) 



Its contribution to the amplitude b sj is given by 6 graphs whose Feynman diagrams are 
similar to those shown in Fig. |^. The 3 diagrams, where the gluon is absorbed by the s-quark, 
read 



{penguin a) 



1 

36 



6 , rrii, 

- + 19 - 12 log — + Gni 

e /i 



a. 



TT 



CFQd{si\Ol}r\b)tr, 



(A.30) 
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while the diagrams where the gluon is absorbed by the 6-quark is given by 



{penguin b) 



— h 5 — 6 log — 

e 



TT 



CFQd{s^\OTR\h)tree 



(A.31) 



Adding the contributions of the two-loop diagrams in eqs. ( |A.22| , |A.23| , |A.24| ) and the 
counterterm contributions in eqs. ( |A . 2 8| , [A730| , |A .311 ), we end up with a total linear term of 
the form 



M2[twoloop + countert] 



TT 



Qd 

72 



6 



- 13 



Qu 



(s7|Of'"|6) 



I tree 



(A.32) 



Of course, this result can be made finite by introducing a corresponding counterterm pro- 
portional to O"^'" given in eq. (|A.25| ) which minimally subtracts the 1/e-pole in eq. ( |A.32|) . 

At this point one might ask why these additional operators appearing in this subsection 
do not appear in the literature [jll| where the singularity structure of these graphs have been 
worked out in order to extract the 0{as) anomalous dimension matrix. We also note that 
there are other new operators induced if one looks e.g. at the analogous two-loop corrections to 
6 — > S7 associated to other Four-Fermi operators. For example, the corrections to O3 generate 
a new magnetic type operator, where nic in eq. ( |A.25| ) is replaced by mf,. 

To absorb all the divergencies one clearly has to enlarge the operator basis in this scheme. 
However, we believe that it is correct to ignore these additional operators for the leading 
logarithmic result for h —>■ sj. The reason is the following: First, the new operators do not mix 
into the old ones given in eq. (|1.2|) at 0{as)', therefore, the old operators run in the same way 
with or without including the new operators. Second, in the 0{a^) matching (at fi = mw), 
the new coefficients have to be zero, as one can easily see. Third, the absence of a log(/i) term 
in eq. ( |A.32| ) indicates, that the four-Fermi operators do not induce any running of the new 
operators at the leading logarithmic level. The only way, the new magnetic operators could run 
is by multiplicative renormalization. But as these operators have zero initial value as discussed 
above, also this effect is unimportant: To leading logarithmic precision the Wilson coefficients of 
the new operators are zero at each renormalization scale. Therefore the new operators certainly 
do not change the leading logarithmic physics. Equivalently, one can say that one can throw 
away the terms proportional to rric in the building blocks in eqs. ( [A.7|) and (|A.8| ) when working 
at leading order. 

However, our result in eq. ([A.32|) shows that the matrix element of O2 leads to finite terms 
(which are of next-to-leading order) linear in mc when calculated in the HV scheme. It would 
be very interesting to see, if these linear nic term are cancelled by the next-to-leading Wilson 
coefficients of the new operators. Such a cancellation is expected to occur, of course, because 
in any other scheme which respects the 4-dimensional chirality properties, these terms do not 
appear. 



B Bremsstrahlung corrections 

In order to make the paper self-contained, we give in this appendix the formulae for the process 
b s'jg based on the operators O2, O7 and Og in the NDR scheme. As we have given the 
analogous virtual corrections to 6 — 57 in the limit 0, we also present the Bremsstrahlung 

corrections for this case. 
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We denote the various contributions to the Bremsstrahlung decay width by T\2, Ty^, Fgg, 
Fgy, F28, Fyg; for example, F22 is based on the matrix element squared of O2, while F27 is an 
interference term between the matrix element of O2 and O7, etc. . 

Note that all the interference terms and F22 are infrared (and collinear) finite. The sum of 
these four finite contributions is given by 

"GSn^mb Jps 

T22 = 2mlQlC'M \Kmml - 2{qr)) (B.2) 

= -32ml Q„ C2{fi)C'/fifi) [qr] Re{t^) (B.3) 

T2S = -32mlQuQdC2{lJ.)Cl" {qr)Re{K) (B.4) 

rrs = -128mlQ,C^,^fi^.)Cl^fi^) (qr) < , (B.5) 

where as in the previous sections p, p', q and r denote the four-momenta of the b- and the 
s-quark, the photon and the gluon respectively. The function k is defined as 

4(2G(t) + 1) 2{qr) 

K = ; t= — (B.6) 

t mi 

G{t) = f^^\og[l-ty{l-y)-ie] . (B.7) 

The phase space boundaries (denoted above by PS) are given by 

i^.e[0,^]; (B.8) 

Note, that the function k in eq. (|B.6|) is finite for m^ 0. Therefore also the Bremsstrahlung 
corrections have a finite limit for m^ 0. 

Fyy is singular for or r || p'. As it cancels the corresponding singularity of the 

virtual corrections to O7, we have taken into account the contribution F^^ already in section 3, 
i.e., F77 is contained in the finite quantity in eq. (|3.3| ). 

Finally, Fgg is singular for Q 01 q \\ jf . These singularities can be removed by adding 

the virtual photon corrections to b ^ sg. This finite sum Fgs can easily be obtained from F77 
in eq. ([3.31). It reads 



'G^g.C|^^A,| f ^ - ^ + 41og^ ) . (B.9) 



88 



To summarize, the total inclusive decay width F for B —>■ + 7 is then given by F = 
Yvirt ^ pfcrems ^ y^s, whcre F"^*, F^,^^™^ and Fgg are given in eqs. (^), (^) and (IROD , 
respectively. 
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